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Chapter I 

Introduction 


1.1. Background 

The motivation for the development of flight test trajectory controllers is well doc- 
umented in the literature [1-6]. In the past, several control system design techniques have 
been applied to this problem, viz, the LQR approach [4], the eigenstructure assignment 
technique [5], Kosut’s suboptimal LQR design with guarenteed stability margin and min- 
imum error excitation criterion [6]. In order to apply these approaches, one is required to 
linearize the aircraft model about the desired flight test trajectory. The linearized aircraft 
model is then a combination of a linear perturbation model and open-loop state-control 
histories along the desired flight test maneuver. A closed-loop perturbation controller can 
then be synthesized using one of the several techniques for linear system design. Since the 
linearized model changes as a function of the flight condition, gain scheduling is invariably 
required to obtain satisfactory control system response. Note that the linearized model 
generates only a part of the total control. The second component of control is provided 
by the trim values of the control variables along the desired trajectory. This part of con- 
trol is open-loop. Thus the flight test trajectory controller iB a combination of closed-loop 
linear perturbation control and the open-loop trim control. While such an approach can 
certainly be made to work for arbitrary nonlinear plants, the amount of data that has to 
be stored can sometimes be prohibitive. In such case, one is forced to compromise and 
store the gains and open-loop controls at sparser intervals. As a consequence, the control 
system suffers performance degradation. It is often possible to partially compensate for 
this performance loss by designing relatively parameter insensitive controllers. However, 
the performance degradation due to lack of fidelity in the trim control settings are often 
more difficult to correct. In addition to this, there is another mechanism responsible for 
the degradation of these controllers, viz, the effect of off-nominal conditions. This arises 
due to the inaccuracy of the linearized model when the aircraft is far from the desired 
trajectory. The trim controls as well as the closed-loop controls in this case would be in 
error. It is interesting to note here that this effect can arise even if the gains and open-loop 
settings are stored at close intervals. Several of these issues were addressed in an earlier 
contract with NASA Ames-Dryden Flight Research Facility [6]. 

Recently, an alternate approach for nonlinear control synthesis has emerged [7-13]. 
The development of this technique for flight controls is due to G. Meyer of NASA Ames 
Research Center. Much of the initial theoretical foundations in this area were laid by 
Brockett [14], Hunt, Su and Meyer [11-13]. In this work, the nonlinear dynamic model is 
assumed to be of the form 


x = f{x) + g{x)u 


( 1 . 1 ) 


The simplest case arises when the number of states in the model are equal to the 
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number of controls. In this case, if g(x ) is invertible and x is perfectly known, the above 
system can be transformed to the form 


where, 


x — V 


( 1 . 2 ) 


V = f{x) + g(x)u, (1.3) 

are the pseudo control variables. A linear feedback controller can be designed for the 
system (1.2). This controller will be of the form 


V = kx. (1.4) 

The expressions (1.3) and (1.4) can now be combined to obtain the nonlinear controller of 
the form 


u = g{x) l [kx - f{x)] (1.5) 

It is clear that if f(x) and g{x) are perfectly known, the nonlinear controller (1.5) will 
have all the attributes of the pseudo linear controller (1.4). More often than not, this 
information is not perfect. As a result, one may have to compensate for this error by- 
building-in extra margins in the controller (1.4) through the introduction of additional 
compensators. Note that in this development, only the control variables were transformed. 
In general, whenever the state-control excess is zero, the prelinearizing transformation will 
leave the original system states unaltered, modifying only the control variables. 

In most practical situations, however, the number of controls are less than the 
number of states. E this case, a simple transformation such as the one described above 
is no longer feasible. A prelinearizing transformation may still be synthesized, but the 
states will not remain unaltered under this transformation. Moreover, the prelinearizing 
transformation would involve the partial derivatives of the right-hand-sides of the original 
nonlinear system. This would require one to impose certain other conditions on the func- 
tions f(x) and g(x). Such requirements are naturally handled in the language of Lie group 
germs and Lie algebras, see Guggenheimer [15] for an initial exposition and Hunt, Su and 
Meyer [11] for specific details. 

E the cases where the number of controls exceed the number of states, one has 
freedom in choosing the combination of controls to be employed in a particular situation. 
This a natural setting for system optimization. One such case has been discussed by 
Cicolani, Sridhar and Meyer [16]. 

Attempting to construct prelinearizing transformations in the case where the num- 
ber of control variables are less than the number of states will invariably lead to the com- 
putation of several partial derivatives, which can sometimes make the problem ill-posed. 
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An alternate approach was proposed recently [17] using singular perturbation theory. Air- 
craft forms a natural system for the application of singular perturbation theory since, 
there clearly are slow and feist modes in the system. For example, the body rates evolve 
faster than other state variables. The use of singular perturbation theory in the aircraft 
control problem can eliminate the need for partial derivatives to a certain extent. This 
can improve the problem conditioning significantly. Moreover, the controller requires less 
computations, making it more efficient for on-line implementation. Reference [17] discusses 
the details of the approach. This paper is given in Appendix A. 

Two modes of operation for the flight test trajectory controllers have been identi- 
fied in the literature [2]. The first is that of trajectory tracking with pilot - in - the - loop 
and the second, a completely automatic operation. The attention will be focussed on the 
latter, though the manual flight test trajectory control aspect will be briefly examined. 

1.2. Summary of Results 

The results of the present study are summarized in the following. 

1. A systematic approach for the development of flight control systems using sin- 
gular perturbation theory and the theory of prelinearizing transformation was 
developed. A paper containing this work has been presented at the Guidance 
and Control Conference at Snow Mass, Co. which is appended at the end of this 
report. 

2. The inclusion of the Command Augmentation System (CAS) in the aircraft model 
modifies the flight test trajectory control problem considerably. Specifically, the 
boundary layer corrections required in the original singular perturbation approach 
are no longer required. This simplifies the control system development consider- 
ably. 

3. The performance of the synthesized controllers are evaluated along the six required 
flight test maneuvers using a complete simulation of the CAS and the aircraft. 

In the next project phase, these controllers will be evaluated in a manned simula- 
tion before attempting a flight test. 

1.3. Report Organization 

A description of the flight test trajectory control system is given in Chapter II. 
This chapter also gives the details of the nonlinear controller synthesis. In Chapter III, the 
controller performance for the six required flight test maneuvers are given. Conclusions 
and future work are outlined in Chapter IV. Appendix A presents a paper describing 
the application of singular perturbation theory for synthesizing nonlinear controllers for 
aircraft. 
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Chapter II 

Flight Test Trajectory Controller Synthesis 


2.1. Aircraft Model 

The six-degrees of freedom model for aircraft flight over a flat nonrotating earth 
is given by : 


V — \—Dco80 4- Ssinfi + {k\cosctcos0 4- kzsinacos0)T — mg(sinBcosotcos/3 
— cosBsin(f>sin0 — cos0cos<f>8inacos/3)]/m 


( 2 . 1 ) 


— [ — L + (& 3 Cosa — k\sina)T 4- mg(cos6cos<j>cosa 4- sinQ since)] /mV cos /3 
4?“ tanfi(pcosa 4- rstna) 


( 2 . 2 ) 


0 = \Dsin0 4- Scosfi — {k\cosasin0 4- kzsinasin0)T 4- mg(sinBcosasin0 
+ cosBsin<f>cos0 — cos6cos<j>sinasin0)]/Vm 4- psina — rcosot 


h — V (cosficosastnO — sin0sin<f>cos6 — cos0sinacos<f>cos0) 
4> = p 4- qsin<f>tan6 4- rcos<f>tanB 
B = qcos<f> — rsin<f> 

p = \PI 1 4* RIz + + C2l\Ba + *3^1 &r 

+ pq{Ix,h ~ D,h) - qr{D x h + U)]// 


(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 


(2.7) 
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q = [Ql 4 + fihh + hUh + hUSr + p 2 h,h ~ prDyh + r 2 I gt I 4 )/I ( 2 . 8 ) 


f = [P /3 + RIq + gi I 6 6 e + g^IeSa + gzI&B r + pq{h,Iz — D,h) 
-qr(D x I 3 + I s ,I 6 )]/I 


( 2 . 9 ) 


With 


h = Iyh 


— lylzM 


h = hh ~ h. 


h = hly 

D x =I t ~ Iy 

Dy = I x h 
I = I x Iyh - IylL 


hy — Iy, = 0 for the Aircraft under consideration. 

The state variables in this model are V the airspeed, h the altitude, a the angle 
of attack, f) the angle of sideslip, 6 the pitch attitude, <f> the roll attitude, p the roll body 
rate, q the pitch body rate, and r the yaw body rate. The yaw attitude 0, down range 
x and the cross range y are ignorable in the flight test trajectory control problem. The 
control variables in this model are throttle, the elevator deflection S e , the rudder deflection 
6 r and the differential tail 6 a . P, Q, and R are the total aerodynamic and thrust moments 
about the roll, pitch, and yaw body axes, not including the moments due to control surface 
deflections. The variables ei, e2» < 3 , / 1 , /a, /s, gi, g-i, gz are the control surface influence 
coefficients for roll, pitch and yaw axes, respectively. Though the aircraft is equipped with 
ailerons, flaps and speed brake, these are accessible only through manual control. The 
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aircraft has a command augmentation system (CAS), which has inputs from the joystick 
and rudder pedals, and also from the autopilot. For automatic flight test maneuvering, it 
will be assumed that the pilot inputs are zero and that the autopilot input port will be 
connected to the flight test trajectory controller. In all that follows, for the sake of brevity, 
the flight test trajectory controller will be termed the Maneuver Autopilot (MAP). For a 
more detailed nomenclature, see the paper given in Appendix A. 

In addition to this, the aircraft engine dynamics of the form 


V = T ( rj c - 17) (2.10) 

will be included in the analysis. In equation (2.10), tj is the actual throttle setting while 
t}e is the commanded throttle setting, r is the engine time constant. In reality, the engine 
dynamics is much more complex than that given by the expression (2.10). But inclusion of 
these in the analysis will complicate the development and may not lead to simple control 
laws. Earlier work [2] has shown that this model is adequate in the flight test trajectory 
controller development. 

If the aircraft did not have a Command Augmentation system, the flight test 
trajectory controller may be designed by splitting the dynamics into slow and fast modes 
and applying the methods of singular perturbation theory and the theory of prelinearizing 
transformations as in [17]. 

On the otherhand, if a command augmentation system is integrated with the 
airframe, the maneuver autopilot should work through this system in order to access the 
aircraft control surfaces. In this case, the formal application of the singular perturbation 
theory is no longer feasible due to the noninvertible CAS dynamics. However, if we assume 
that the CAS dynamics is stable and sufficiently fast, one needs to construct only the slow 
controller, leaving the task of boundary layer corrections to the CAS. This approach will 
be followed in the present work. 


2.2. Command Augmentation System 

The CAS forms the interface between the pilot, the autopilot, and the aircraft 
control surfaces. Since the present effort is in generating a completely automatic maneuver 
autopilot, the MAP has to effect trajectory tracking through the autopilot inputs in the 
CAS. 


The autopilot inputs of the command augmentation system consists of three con- 
trol channels, the input commands of which are a n the normal acceleration in g' s, p the roll 
rate and /? the angle of sideslip. It is important to note that these inputs were not supplied 
by the aircraft manufacturer, but were identified after several simulations on the aircraft 
and CAS system. We assume here that the CAS system is capable of controlling these 
three quantities with a sufficient degree of accuracy and speed. Thus, if a MAP system 
generates these input quantities, we leave the burden of maintaining these at the desired 
value to the CAS system. It is to reasonable to assume that the CAS has stable, fast 
dynamics because, otherwise, the aircraft would have been difficult to difficult to fly in the 
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first place. If the resulting system is found to have an unacceptable performance, it may 
become necessary to provide additional feedback compensation for the CAS. Fortunately, 
this case did not arise in the present research. 

The CAS generates the control surface deflections based on the commanded values 
of the normal acceleration in j’s, the angle of sideslip, and the roll rate. The feedback 
quantities for the CAS pitch channel are roll and pitch angular acceleration, pitch body 
rate, normal acceleration and angle of attack, and a stall limit. This channel has gains 
scheduled as functions of Mach number, static and total pressure. In the roll CAS channel, 
the autopilot provides the roll rate command while the feedback quantity is the actual 
aircraft roll rate. The yaw CAS channel command is the angle of sideslip 0 while the 
feedback quantities are the aircraft lateral acceleration, yaw and roll angular acceleration, 
and yaw and roll rates. This channel has a gain scheduled as a function of angle of attack. 
All the control surface actuators have position and rate limits in addition to freeplay and 
deadband nonlinearities. 


2.3. Controller Synthesis 

With the foregoing, in the following, the nonlinear controllers for this system will 
be developed. Since the CAS tracks the commanded values of normal acceleration a n , 
in 0’s, roll rate p and the angle of sideslip /?, the MAP will be designed to generate the 
required commands in these quantities in order to track the commanded values of altitude, 
airspeed, angle of attack and roll attitude. 

The nonlinear controllers will be developed for each of these quantities separately. 
The MAP for a particular maneuver may be synthesized by employing the appropriate 
combinations of these controllers. 


2.3.1. Controller for the Airspeed V 

The state equations for these quantities may next be used to develop this controller. 
Firstly, we shall consider the airspeed loop. The control variable in this case is the throttle, 
the other quantities on the right hand side of the state equation (2.1) being computed from 
the measurement of states. The prelinearizing transformation in this case is particularly 
simple, Viz, 


With, 




( 2 . 11 ) 


U\ —[— Dcos0 + SsinfJ + ( k\cosaeos0 + kssinacos0)T — mg[sinOcosacos0 
— cos6siruj>sinacos0)]/m 


( 2 . 12 ) 
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The aircraft drag D and the sideforce S can be computed using a simple nonlinear model 
or interpolated from stored tables. The thrust component coefficients ki and k^ are known 
constants. Since other terms in (2.12) can be computed from feedback values of the states 0 , 
<t>> /?, if the pseudo control U\ is known, this expression may be inverted to obtain thrust 

T. The thrust may then be reverse interpolated to yield the throttle setting using a table 
lookup. The transformed linear dynamic model for airspeed is given by the expression 
(2.11). A simple proportional controller earn stabilize the system. However, Bince the 
controller has to track ramp airspeed command and because the nonlinear functions on 
the right hand side of expression (2.12) cannot be computed exactly, an integral feedback 
will be used in this channel. With this, the pseudo control U\ will be of the form 

U\ = G\ey + Gj / eydt (2. 13) 

Jo 


Where 


«v = V e - V (2.14) 

V e is the desired value of airspeed, G i and G 2 are the feedback gains to be determined 
based on the maximum available thrust, and the desired speed of response. These gains 
may be determined without difficulty since the closed-loop transfer function of the airspeed 
control loop in the pseudo control variable U\ is given by 


V (s) G\8 -f- G 2 

V c (fi) S* + GrjS + Crj 


(2.15) 


The gains G\ and G 2 can be selected to provide adequate damping and speed of 
response. Note, however, that the dynamics given in (2.15) does not explicitly include the 
aircraft thrust limits. It is the designer’s responsibility to pick the combination of gains 
that would lead to a satisfactory performance without violating the thrust limit. 

Using equations (2.13) and (2.12), the actual control can be computed as 


mil 1 + Dcos/3 — Ssinfi + mg(sin$cosacos(J — co8 0sin<f>sinacos0) 

(kicosacosfi + k$stnacosl3) 

The thrust emerging from the expression (2.16) can next be used to determine the desired 
throttle setting r) c from a table. The thrust-throttle relationship for aircraft are often 
highly nonlinear, requiring an iterative algorithm to compute the throttle from the given 
thrust. In the present case, a modified linear interpolation method was used, see page 
10 of Ref. 18 for deatils. In all the maneuvers studied here, this thrust-throttle iteration 
converged in four steps for an absolute thrust tolerance of 1 x 10 -13 . In practice, a higher 
value of tolerance can be used. The computed throttle setting must satisfy 


= T (2.16) 
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Vmin ^ 1}c ^ Vmaz 


( 2 . 17 ) 


If the calculated throttle Betting is larger than the maximum or less than the minimum, 
it may be fixed at the appropriate bound. Note that no explicit compensation for the 
engine time lag has been incorporated. If required, a lead-lag engine prefilter may be 
used to partially offset the engine lag effects. However, this was not found necessary for 
satisfactory operation in the present case. For the aircraft engine considered here, there 
is a small core thrust saturation region approximately between 83° and 98° of throttle 
setting. Since the slope of the throttle-thrust curve is ill defined in this region, a hysterisis 
safe guard logic was introduced. According to this logic, when the thrust is decreasing 
from the maximum afterburner value to that corresponding to less than 98° throttle, the 
throttle is held at this value, till the required thrust falls below that corresponding to 
83° throttle. Reverse of this logic applies whenever the required thrust approaches the 
saturation value from the core thrust direction. 

In some flight test trajectories such as the constant thrust windup turn, the throt- 
tle would have to be fixed while regulating the airspeed. This can be accomplished by 
introducing an altitude rate in the system. Since, 


V = [—Dcosfl + Ssinfl + {k\cosacos0 + k^sinacos0)T — m(y/V)A]/m (2.18) 

In this case, the right hand side of the expression (2.18) would be the pseudo control U\, 
and the actual control would be the altitude rate h. The right hand side of expression 
(2.18) may then be equated to expression (2.13) to compute the altitude rate required to 
maintain the airspeed at the desired value. Therefore, 


(V/mp )[— Dcosfi + Ssinfl -f (kicosacosfi + A* sinacos0)T — mU\ — h (2-19) 

The task of generating this desired altitude rate would be that of the altitude 
control channel. In this situation, it is clear that one would not be able to track an 
arbitrary altitude history, unless it happens to be the same as that emerging from the 
controller. Note that this controller will perform satisfactorily only if the altitude control 
loop can produce the desired altitude rate sufficiently fast. Thus care should be taken 
to ensure that the airspeed control loop remains slower than the altitude loop in this 
maneuver. 

2.3.2. Controller for Altitude h 

Just as in the airspeed channel, two versions of this controller will be required to 
execute the flight test trajectories under consideration. This arises because in certain flight 
test maneuvers such as the excess thrust windup turn flight test trajectory, one would be 
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required to track an angle of attack history while maintaining the altitude at the desired 
value. 


In order to develop the altitude controller, we shall use an alternate form of the 
state equation than that in (2.4). It can be verified that 

h — a x sinO — a v sin<f>cosO — a z cos<f>cos6 (2.20) 

In the above expression o*, a„, and a, are the acceleration components along the X,Y, 
and Z body axes. Now, the normal acceleration o n is related to a* as 


fl* = -On + ( g / go)cos0cos4> (2.21) 

With this, the equation (2.20) becomes, 

m ..... " s - 

h — a x sinO - a v sin<f>cos6 + a n cos<f>cosO - (g/go)cos 2 0cos 2 <t> (2.22) 

Equation (2.22) is in the form suitable for the development of a nonlinear control law for 
altitude. The expression (2.22) may be written as 


h = U 2 


(2.23) 


With, 


U 2 = a x sin6 — a v 8\n<j>cos6 + a n cos<f>cos$ — (g /go)cos 2 Ocos 2 
A linear control law of the form 


(2.24) 


r f 

U 2 = + G^e'h + Gs / e^dt 

Jo 


(2.25) 


Where 


eh = (h e - h) (2.26) 

may be set up to obtain desired time response for the transformed system (2.23). The 
calculated U 2 may then be used to compute the normal acceleration command to the CAS 

as ■ -■ ' " - ■ 


U 2 — a x sin6 + a y sin<f>cos6 + (g/go)cos 2 0cos 2 <f> 
cos<j>cos9 


(2.27) 
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As in the airspeed control loop, the closed-loop transfer function for the altitude 
control loop in the pseudo control variable can be obtained as 

h(g) _ G3S + G4S 2 + G$ 2g) 

h c (s) a® 4- G$s ■+- -+- G$ 

The choice of the gains G 3 , G 4 , and G 5 are a little more difficult in this case, 
because the expression (2.28) is a third order transfer function. However, this is straight 
forward. 

As mentioned earlier, in certain flight test trajectories, one would be required to 
control the angle of attack and the altitude simultaneously. In this situation, the normal 
acceleration is nolonger available as the control variable to maintain the altitude at the 
desired level. The above approach will then have to be modified. 

Whenever the simultaneous control of angle of attack and angle of attack are 
desired, one can assign the roll attitude <f> as the control variable in the altitude channel. 
If this is done, the roll attitude required to control the altitude can be computed as 


4>c = COS 


-\ 


a x sinO — U% 


■cosO 




l+>v 


+ < l > 0 


(2.29) 


Where, 


^0 — tan 1 



(2.30) 


The assumption here is that the dynamics of the roll attitude control loop is much faster 
than that of the altitude control loop, a reasonable assumption for high performance fighter 
aircraft. This commanded roll attitude <f> e must be tracked by the roll attitude control loop. 
Note that the commanded roll attitude emerging from the expression (2.29) is devoid of 
any sign. This indicates that the altitude control may be achieved through a roll attitude 
yielding either a right handed or left handed turn. In the present work, it will be assumed 
that positive values of <f> c will be used, leading always to a right-hand turn. 

Before closing this section it is perhaps worthwhile to note that one should be 
careful in selecting too high altitude loop gains. Clearly, the aircraft has a maximum nor- 
mal acceleration limit beyond which the control system performance cannot be predicted. 
Moreover, when the control law (2.29) is being used, trying to speed up the altitude control 
channel beyond a certain point will result in large changes in the commanded roll attitude 
<t> c which the roll attitude controller may not be able to track. This can lead to control 
system instability. 


2.3.3. Controller for Angle of Attack a 

A nonlinear controller for the angle of attack a can be set up in a manner analogous 
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to that of the airspeed control loop. As before, we define the right-hand-side of expression 
(2.2) as another pseudo control variable U 3 . Thus, 


With 


d = U 3 


(2.31) 


Uz — [— L + (fc 3 cosa — kisina)T -f- mg{cos6cos4>cosa + sinBsina)]/(mVcosP) 


+9 - tan(3(pcosa + rsina) 
A linear control law of the form 


Us — Geto + G-j 



e a dt 


Where, 


(2.32) 


(2.33) 


e a = a c -a (2.34) 

may be designed to obtain the desired time response from the prelinearized dynamic system 
(2.31). In expression (2.34), a c is the commanded value of angle of attack. The gains G& 
and G 7 can be chosen based on the desired natural frequency and damping ratio. The 
closed-loop transfer function of the angle of attack loop in the pseudo controls is given by 


o(i) Ggs + Gj 

ac c (s) s 2 + Ggs + Gf 


(2.35) 


The actual control variable in the angle of attack channel is the normal accelera- 
tion. This variable can be obtained by using the expression (2.32) in conjunction with the 
expression for normal acceleration. Since the normal acceleration is given by 


= [— kzT -I- Dsina + Lcosa — mgcos6cos<f>]/mg, (2.36) 

the expression (2.32) and (2.35) may be used to obtain the normal acceleration required 
to track a given angle of attack history as : 
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= [ y ”fg £ g g] | \DUma - kiTsina 


+mg{co30eos^co3a + ain05tn<l)]/(mV cosp) + g — tan0[peo3a + miw) — Vz j — cosdcos# 

(2.37) 

This completes the development of the angle of attack controller. 

2.3.4. Controller for the Roll Attitude <f> 

The controller synthesis for the roll attitude <)> follows the same steps as those 
described for airspeed and angle of attack channels. We first assign a pseudo control 
variable C /4 such that 


4> = U 4 


(2.38) 


With 


U4 = p + qsin<f>tan8 + rcos<f>tan8 (2.39) 

We next design a linear controller of the form 


U4 = G&Cj, + Gg 



(2.40) 


Where, 


e* = tc ~ <t> ( 2 - 41 ) 

The variable <f> c is the commanded value of roll attitude in expression (2.41). The closed- 
loop transfer function for the roll attitude control loop in the pseudo control variable is 
given by 


^(*) _ Gs* + Gg, ( 2 . 42 ) 

^ c (s) S 2 + GgS + Gg 

It is clear that the gains Gg and Gg can be chosen to meet the desired time response spec- 
ifications. The pseudo control variable obtained from (2.40) may be used in the expression 
(2.39) to obtain the actual control variable p. Thus, 

p = U4 — qsintptanO — rcos<f>tanB (2-43) 
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It is clear that as long as the fedback state variables q, r, <j>, and 6 are close to their 
actual values, the nonlinear controller will have almost the same performance as that of 
the pseudo control loop. 

This completes the closed-loop controller development. These control laws may 
be further simplified because all the maneuvers are to be carriedout with (3 — 0. We note 
here that in order to implement these control laws, one needs to be able to compute the 
aircraft drag and the throttle setting from the calculated thrust. There are two ways to 
approach this problem. The first one consists of storing these as tables and the second, 
to construct simple approximate models for these quantities. The latter approach will be 
used in the present work. The state variables required are V, h, a, 6, <f>, q, r, and the 
three accelerations a*, a y , and a M . In the next Chapter, the controllers developed here 

will be combined together to form the maneuver autopilot to execute particular flight test 
trajectories. 

In the following section, we shall briefly indicate how the synthesized closed-loop 
control laws may be modified for operation with pilot-in-the-Ioop. 

2.S.5. Maneuver Control with Pilot-in-the-loop 

In the foregoing, the technique for the synthesis of nonlinear control laws for flight 
test trajectory control was discussed. In order to enable manual control, these control laws 
may be modified in several different ways. One approach would be assume that the errors 
ev ' e *> e 4> remain constant over a time interval to predict the pseudo controls U\, 

U 2 , U 3 , and U 4 . These control variables may be displayed to the pilot through devices 
such as the one described in [2]. This approach compensates only the observation and 
neuromuscular delay in the human pilot dynamics [4]. 

An alternate approach would be to include the pilot transfer function in the pseudo 
control loops and select the gains to achieve the desired time response. Mere gain stabi- 
lization may not always be feasible in this case. If the resulting control loop displays 
unacceptable behavior, appropriate dynamic compensation networks will then have to be 
introduced. In any case, the present methodolgy give a systematic approach to handle the 
manual control problem also. 



Chapter III 

Contoller Evaluation 


3.1. Simulation 

In order to evaluate the performance of the nonlinear flight test trajectory con- 
troller, a Fortran implementation of this controller is mechanized on a complete CAS -(-Aircraft 
model developed at NASA Ames-Dryden Flight Research Facility. This section describes 
the details of controller implementation and the results of closed loop simulations with the 
nonlinear controller for the six required flight test manuevers. 

Specifically, the following maneuvers are simulated. 

- Level Acceleration Trajectory 

- Pushover Pullup Trajectory 

- Zoom and Pushover Trajectory 

- Excess Thrust Windup Turn Trajectory 

- Constant Thrust Windup Turn Trajectory 

- Constant Dynamic Pressure, Constant Load Factor Trajectory 
For a detailed description of these maneuvers, see Ref. 17. 

The simulations are mechanized on the SYSTEM-BUILD software of Integrated 
Systems Inc., some details of which are given in [19,20]. This proprietory software is block 
diagram oriented and is suitable for efficient generation of nonlinear, multi-rate simulations. 
This software incorporates several commonly encountered dynamic subsystems in a library. 
The cases where the dynamics is not representable with any of the available blocks in the 
library, Fortran blocks may be attached. This is often the situation in aircraft simulations 
wherein the aerodynamics and engine thrust are given as nonlinear functions of several 
state variables. In the next section we give a brief description of the aircraft and CAS 
simulation as it is mechanized on SYSTEM-BUILD. 

3.2 SYSTEMJBUILD simulation of Aircraft and Controller 

Airframe, Engine and command augmentation system models provided by NASA 
Ames-Dryden Flight Research Facility are linked with the Fortran code for the nonlinear 
flight test trajectory controller in the SYSTEM-BUILD environment. A block diagram 
of this implementation is given in Figure 3.2.1. In this figure, the CAS,AROMDL and 
ENGINE super blocks in the closed loop simulation were supplied NASA Ames-Dryden 
Flight Research Facility, while the CONT super block was developed and linked to the 
SYSTEM-BUILD model of the aircraft and CAS. The listing Fortran code implemented 
in the block CONT is given in Appendix B. 


15 



Figure 3.2.1. SYSTEM-BUILD Block Diagram of the Simulation. 
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The command augmentation system is implemented as as a discrete Fortran block 
in the simulation. The sampling interval of 0.02 seconds has been fixed from previous sim- 
ulations of the aircraft and CAS system. The nonlinear aircraft model is also implemented 
as Fortran block, but it is a continuous system. 

Real time simulation and testing of the nonlinear flight test trajectory controller 
will require that the output of the controller be interfaced to the aircraft through the CAS 
system. This in turn requires that the output of the controller be the normal acceleration, 
angle of sideslip and the roll body rate rather than control surface deflections. Since the 
CAS nonlinear dynamics are not invertible, one has to either assume that the CAS is stable 
and sufficiently fast when compared with the MAP such that it tracks the input commands, 
or build-in a feedback compensation to ensure that it tracks the input commands. The 
present experience indicates that the former assumption is valid. 

r ^ me scale separation is a key feature of the nonlinear controller, and here, we 
briefly explain how this separation was implemented in the computer simulation of the 
closed - loop aircraft model. Given that the CAS model has been coded to run in the 
SYSTEM-BUILD simulation as a discrete time block with a fixed step size of 0.02 seconds, 
a second order Runge-Kutta algorithm was used for all simulations with a step size of 
0.02 seconds. This time scaling is consistent with the dynamics of the continuous time 
airframe and engine models but much too fast for the nonlinear controller. In order to 
slow the command variables from the controller so that response times of the aircraft and 
CAS states were shorter than the dynamics of the controller commands, a flag internal to 
the controller code updates the command outputs every V integration steps. Response 
testing at straight and level flight conditions with »' = 5 gave good results. Thus, the 
foilwing simulations were all completed with the nonlinear controller running at 10 Hz. 


3.3. Flight Test Trajectory Control Simulation 

Six required flight test trajectories were simulated to evaluate the performance of 
the nonlinear controller. One of the major problems encountered during the simulation 
was that of finding a consistant set of initial conditions that would trim the aircraft at the 
desired altitude and airspeed. Though the NASA-DFRF linearization program produces 
consistant set of initial conditions, it was found that when these were introduced in the 
simulations, there was a significant starting transient. For now, this transient is attributed 
to the unknown initial conditions in the CAS , since this system contains several dynamic 
loops. Moreover, some repeatability problems were encountered in the initial phases which 
were traced to two undefined logical variables in the CAS model. As mentioned earlier, 
the CAS input variables were not clearly known at the outset, but had to be identified 
after several simulations. 

In the following, the performance of the nonlinear controller along each of the six 
flight test maneuvers will be given separately. 

3.3.1. Level Acceleration Flight Test Trajectory 

The requirements in this maneuver are : track a ramp airspeed history while 
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maintaining constant altitude with wings level. To achieve this objective, the nonlinear 
controller for the airspeed given by the expressions (2. 13), (2.14), and (2.16) is used in 
conjunction with the altitude controller expressions (2.25), (2.26) and (2.27), and the roll 
attitude controller expressions (2.40), (2.41) and (2.43). The commanded roll attitude <f> c 
is zero in this case. 

A level acceleration flight test trajectory starting at 30000 feet and 0.9 Mach is 
given in Figures 3.1 through 3.4. In this maneuver, the aircraft was required to accelerate 
from 0.9 Mach to 1.2 Mach in about 40 seconds while maintaining the altitude within ± 
50 feet. This represents a particularly stringent maneuver since the aircraft has to pass 
through the transonic region while maintaining the altitude. From Figure 3.1, it can be 
seen that the aircraft tracked the Mach number history with very little error, except in 
the vicintiy of Mach 1. In this figure, the dotted line represents the commanded Mach 
number, while the solid line shows the tracked history. If a linear perturbation controller 
had to perform this task with an equivalent accuracy, it would have required atleast three 
sets of gains, the first set in the subsonic regime, the second in the transonic regime and 
the third one in the supersonic regime. The altitude history given in Figure 3.2. shows 
that the altitude was maintained within ±15 feet, well within the require accuracy. The 
elevator history given in Figure 3.3 shows about 1 HZ oscillation in the transonic region 
getting rapidly damped out as the aircraft commenses its supersonic flight. However, the 
elevator deflection is well within the saturation limits, the maximum deflection being about 
^®°* ■^ roni Figure 3.4, it may be observed that the throttle saturates during the transonic 
region at its maximum afterburner value and continues to stay there for about 5 seconds. 
Subsequently, oscillations can be seen which rapidly damp out during the supersonic flight. 
Overall response of the control system is very good. Considering that this performance is 
achieved with a single set of gains, the performance is indeed remarkable. 


3.3.2. Pusfiover/Pullup Trajectory 

The requirement here is to track a saw tooth wave form in angle of attack while 
maintaining the constant Mach number. The maneuver autopilot for this trajectory con--~ 
sists of the airspeed controller given by the expressions (2.13), (2.14) and (2.16), the angle 
of attack controller given by the expressions (2.33), (2.34) and (2.37), and the expressions 
(2.40), (2.41) (2.43) for the roll attitude controller . 

Figures 3.5 though 3.9 give the behavior of the controller along a typical pushover /pullup 
flight test trajectory. The commanded angle of attack history is shown in dotted line in 
this figure while the control system response is Bhown in solid line. The small amplitude 
oscillations in angle of attack at the begining of this maneuver is due to the initial condi- 
tions rather than due to the command. After the initial oscillations, the angle of attack 
history tracking is very good. The Mach number is maintained within 0.004 throughout 
the maneuver. The altitude is of no particular concern here, since it is an uncontrolled 
variable. However, for the sake of completeness, the resulting altitude is given in Figure 
3.9. The controller performance is again very good. 


18 



1.25 


1.2 


1.15 

M 

A 1.1 
C 

H 1.05 
1 

.95 


N 

0 


.9 



Figure 3.1. Mach Number vs. 

Dotted line: 1 
Solid line: 



TIME , S 



19 






TIME , S 


Figure 3.3. Elevator Deflection vs. Time along the Level Acceleration Trajectory. 
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Figure 3.5. Angle of attack vs. Time along the Pushover/Pullup Maneuver. 

Dotted Line: Command Angle of Attack. 

Solid Line: Actual Angle of Attack. 
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Figure 3.8. Elevator vs. Time along the Pushover/Pullup Maneuver. 
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Figure 3.9. Throttle Setting vs. Time along the Pushover/Pullup Maneuver. 
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3.3.3. Zoom and Pushover Flight Test Trajectory 

The zoom and pushover flight test trajectory requires flight along a parabolic flight 
path with a set of Mach/angle of attack/ altitude condition to be met at the apex of the 
parabola. The flight along the parabolic path occurs with constant thrust. A systematic 
approach for generating commands for Buch a flight test trajectory was discussed in the 
earlier contract phase [6]. In the present work, however, it was found that simply by 
stretching the angle of attack history in the pushover/pullup maneuver, one could obtain 
the zoom and pushover trajectory. The maneuver autopilot for this trajectory consists of 
the expressions for the airspeed controller (2.13), (2.14) and (2.16), the angle of attack 
controller expressions (2.33), (2.34) and (2.37), the roll attitude controller expressions 
(2.40), (2.41) and (2.43) in the initial transient phase. During the zoom and pushover 
phase, the airspeed control loop is disabled since this maneuver is required to be flown at 
constant throttle. 

The desired conditions at the apex of the parabola in the present case was 30000 
feet altitude, 0.7 Mach and 3° angle of attack. The actual values at the apex were, 30019.26 
feet, 0.691 Mach and 3.17° angle of attack. The actual errors are well within the required 
specifications. 

The trajectory variables for this maneuver are presented in Figures 3.10 through 
3*14. The throttle was fixed at 20% during the zoom and pushover maneuver. The 
controller performance is found to meet the desired accuracy specifications. 


3.3.4.Excess Thrust Windup Turn 

The requirement here is to track a linear angle of attack history while maintain- 
ing the airspeed and altitude constant. This maneuver is highly coupled and potentially 
unstable because the aircraft can have bank angles very close to 90°. 

The nonlinear controller uses the normal acceleration to track the desired angle of 
attack history, the throttle to maintain the airspeed and the roll attitude to maintain alti- 
tude. Thus, the maneuver autopilot consists of the airspeed controller expressions (2.13), 
(2.14), (2.16), the angle of attack controller expressions (2.33), (2.34), (2.37), the alti- 
tude controller equations (2.25), (2.26), (2.29), (2.30) in conjunction with the roll attitude 
controller equations (2.40), (2.41) and (2.43). = ; 

Alternately, the maneuver modeling scheme described in [6] may be used to gen- 
erate a consistent set of angle of attack and roll attitude commands which may be tracked 
using the V , a, <f> controllers to obtain the desired excess thrust windup turn trajectory. 
This approach is employed here in the interest of simplifying the controller. Note that with 
this approach, the altitude is not controlled, but is allowed to vary as it may throughout 
the maneuver. 

The performance of the nonlinear controller along this maneuver is illustrated in 
Figures 3.15 through 3.26. There appears to be hang-off error in the angle of attack channel 
which increases with an increase in the roll attitude. Since there is an integral feedback 
in this channel, such a hang-off error can arise only if the CAS input in the pitch channel 
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Figure 3.11. Angle of Attack vs. Time along the Zoom and Pushover Maneuver. 

Dotted Line: Commanded Angle of Attack 
Solid Line: Actual Angle of Attack. 
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Figure 3.12. Mach number vs. 



TIME , S 


ime along the Zoom and Pushover Maneuver. 
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Figure 3.13. Elevator Deflection vs. Time along the Zoom and Pushover Maneuver. 
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Figure 3.14. Throttle Setting vs. Time along the Zoom and Pushover Maneuver. 
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Figure 3 . 15 . Angle of Attack vs. Time along the Excess Thrust Windup Turn Flight Test 
Trajectory. 


Dotted Line: Command Angle of Attack. 
Solid Line: Actual Angle of Attack. 
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3.16. Mach Number vs. Time along the Excess Thrust Windup Turn Flight Test Tra- 
jectory. 
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Figure 3.17 Roll Attidute vs. Time along the Excess Thrust Windup Turn Flight Test Trajec- 


tory. 


Dotted Line: Commanded Roll Attitude. 


Solid Line: Actual Roll Attitude. 
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Figure 3.19. Angle of Sideslip vs. Time along the Excess Thrust Windup Turn Flight Test 
Trajectory. 
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3.20. Roll Body Rate (Deg/s) vs. Time along the Excess Thrust Windup Turn Flight 
Test Trajectory. 
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Figure 3.21. Pitch Body Rate (Deg/s) vs. Time along the Excess Thrust Windup Turn Flight 
Test Trajectory. 


39 



0 2 4 6 8 10 12 14 16 18 20 

TIME, S 


Figure 3.22. Yaw Body Rate (Deg/s) vs. Time along the Excess Thrust Windup Turn Flight 
Test Trajectory. 


40 






0 2 4 6 8 10 12 14 16 18 20 

TIME, S 

Figure 3 . 24 . Elevator Deflection vs. Time along the Excess Thrust Windup Turn Flight Test 
Trajectory. 
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Figure 3.25. Differential Tail Deflection vs. Time along the Excess Thrust Windup Turn Flight 
Test Trajectory. 
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Figure 3 . 26 . Rudder Deflection vs. Time along the Excess Thrust Windup Turn Flight Test 
Trajectory. 
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was a component of the normal acceleration rather than the pure normal acceleration. It 
is possible that this hang-off error arises due to angle of attack limit also. In any case, 
the angle of attack error is less than a degree throughout the maneuver. The important 
requirement in this maneuver, however, is that the angle of attack history be linear in 
time, which is certainly satisfied by the present nonlinear controller. The Mach number 
error is within 0.0005 throughout the maneuver. The roll attitude tracking is excellent, 
being less than 0.2 degrees throughout the maneuver. The altitude error throughout the 
maneuver is within ±50 feet, which satisfies the specifications. It should be possible to 
decrease this error by using finer trim data in the maneuver modeling program. 


3.3.5. Constant Thrust Windup Turn 

In this maneuver, it is desired to track a linear angle of attack history with fixed 
throttle setting while maintaining constant Mach number. Since the throttle is fixed, 
the Mach number can be maintained only through an altitude rate. Moreover, since the 
normal acceleration channel is used to track the angle of attack history, the required 
altitude rate will have to be achieved through an appropriate roll attitude. Thus the 
maneuver autopilot consist of the airspeed controller (2.13), (2.14), (2.19) in conjuction 
with the altitude controller (2.25), (2.26), (2.29), (2.30) and the roll attitude controller 
(2.40), (2.41), and (2.43); and the angle of attack controller (2.33), (2.34) and (2.37). Note 
that this implementation relies on the fact that the airspeed control loop is driven slower 
than the altitude control loop which in turn is assumed to be slower than the roll attitude 
control loop. 

Just as in the excess thrust windup turn flight test trajectory, in order to simplify 
the controller implementation, the maneuver modeling described in [16] is used to generate 
a consistent set of angle of attack and roll attitude commands to yield the desired constant 
thrust windup turn trajectory. 

The performance of the controller along the constant thrust windup turn trajectory 
is given in Figures 3.27 through 3.38. Just as in the excess thrust windup turn case, there 
is a hang-off error in the angle of attack history, which appears to depend on the roll 
attitude. The tracking error is always less than about one degree throughout. The Mach 
number is maintained within 0.015 during the maneuver while the altitude is within 120 
feet of the required value. Throughout this maneuver, the throttle is maintained at 105 
degrees. 


3.3.6. Constant Dynamic Pressure, Constant Load Factor Trajectory 

Along this flight test trajectory, the load factor and the dynamic pressure should 
be constant, while maintaining a desired Mach rate. The load factor can be maintained 
by fixing the value of normal acceleration while the desired Mach rate can be sustained 
using the throttle. The dynamic pressure is maintained by an appropriate altitude rate 
depending on the desired Mach rate and ambient density. The altitude rate is generated 
through the vehicle roll attitude. It is clear that one cannot permit an arbitrary Mach 
rate in the system while maintaining the load factor and dynamic pressure constant since 


45 




TIME , S 


Figure 3.27. Angle of Attack vs. Time along the Constant Thrust Windup Turn Flight Test 
Trajectory. 

Dotted Line: Commanded Angle of Attack. 

Solid Line: Actual Angle of Attack. 
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Figure 3.28. Mach Number vs. Time along the Constant Thrust Windup Turn Flight Test 
Trajectory. 
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Figure 3.29. Roll Attitude vs. Time along the Constant Thrust Windup Turn Flight Test 
Trajectory. 

Dotted Line: Commanded Roll Attitude. 

Solid Line: Actual Roll Attitude. 
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Figure 3.31. Angle of Sideslip vs. Time along the Constant Thrust Windup Turn Flight Test 
Trajectory. 
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Figure 3.32. Roll Body Rate (Deg/s) vb. Time along the Constant Thrust Windup Turn Flight 
Test Trajectory. 
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Figure 3.34. Yaw Body Rate (Deg/S) vs. Time along the Constant Thrust Windup Turn Flight 
Test Trajectory. 
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Figure 3.36. Elevator Deflection vs. Time along the Constant Thrust Windup Turn Flight Test 
Trajectory. 
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Figure 3.37. Differential Tail Deflection va. Time along the Constant Thrust Windup Turn 
Flight Test Trajectory. 
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Figure 3.38. Rudder Deflection vs. Time along the Constant Thrust windup Turn Flight Test 
Trajectory. 
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the atmospheric density gradient is fixed. For a detailed analysis of this maneuver, see [6]. 
The maneuver autopilot for this trajectory consists of the airspeed controller given by the 
expressions (2.13), (2.14), (2.16), and the altitude controller given by (2.25), (2.26), (2.29), 
(2.30) in conjunction with the roll attitude controller (2.40), (2.41) and (2.43). 

Similar to the excess thrust windup turn and constant thrust windup turn flight 
test trajectories, in the present work, this maneuver is implemented using the controllers 
for airspeed, angle of attack and roll attitude for simplicity. 

The controller performance along the constant dynamic pressure, constant load 
factor trajectory is given in Figures 3.39 through 3.51. From Figure 3.41, it can be seen that 
the dynamic pressure has been maintained constant within 0.8% throughout the trajectory. 
The flight test trajectory begem at 10 seconds in the present case, and the load factor is 
maintained within 0.05 of the required value, as can be seen from Figure 3.42. As noted in 
the two earlier maneuvers, there is a hang-off error in the angle of attack channel. During 
the initial transient region, there is a saturation in the elevator and differential tail due 
to the control deflection authority limits. The overall performance of the controller is well 
within the performance specifications. 

S.4. Conclusions 

The closed loop simulation reveals that the initial assumption with regard to the 
speed of the CAS is valid. This system is indeed stable and sufficiently fast. The second 
conclusion that emerges from the present analysis is that the nonlinear Maneuver Autopilot 
is sufficiently robust with respect to modeling inaccuracies since a simpler model than that 
implementd in the simulation could control the aircraft with a high degree of accuracy. 
If the CAS and aircraft models are close to the actual, one would expect a satisfactory 
control system response during the actual flight test also. 
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Figure 3.40. Mach number vs. Time along the Constant Dynamic Pressure /Constant Load 
Factor Flight Test Trajectory. 

Dotted Line: Commanded Mach Number. 

Solid Line: Actual Mach Number. 
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Figure 3.41. Dynamic Pressure vs. Timealonj 
Factor Flight Test Trajectory. 
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Figure 3.42. Load Factor vs. Time along the Constant Dynamic Pressure /Constant Load 
Factor Flight Test Trajectory. 
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Figure 3.43. Angle of Attack vs. Time along the Constant Dynamic Pressure/Constant Load 
Factor Flight Test Trajectory. 


Dotted Line: 
Solid Line: 


Commanded Angle of Attack. 
Actual Angle of Attack. 
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Figure 3.44. Roll Attitude vs. Time along the Constant Dynamic Pressure/Constant Load 
Factor Flight Test Trajectory. 

Dotted line: Comman ded R oll Attitude. 

Solid line: Actual Roll Attitude. 
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Figure 3.46. Pitch Body Rate (Deg/s) vs. Time along the Constant Dynamic Pressure/ Constant 
Load Factor Flight Test Trajectory. 
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Figure 3.49. Elevator Deflection vs. Time along the Constant Dynamic Pressure/Constant 
Load Factor Flight Test Trajectory. 
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Figur 3.50. Differential Tail Deflection vs. Time along the Constant Dynamic Pressure/Constant 
Load Factor Flight Test Trajectory. 
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Figure 3.51. Rudder Defieciton vs. Time along the Constant Dynamic Pressure/Constant Load 
Factor Flight Test Trajectory. 
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Chapter IV 

. Conclusions and Future Work 


This report dealt with the development of nonlinear flight test trajectory con- 
trollers for a high performance fighter aircraft. This work employed the singular pertur- 
bations arguments to reduce the order of the system and the recently developed theory 
of prelinearizing transformations to generate explicit nonlinear controllers. Appropriate 
dynamic compensators were incorporated in this controller to ensure a measure of insensi- 
tivity to modeling inaccuracies. The modifications of these nonlinear controllers to permit 
manual control have been indicated. This would involve the use of a lead-lag filter to 
compensate for the dynamics of the human pilot. 

The closed-loop performance of the nonlinear controller along the six desired flight 
test trajectories were demonstrated through a complete simulation of the aircraft and CAS 
system along with the controller. The controller performance has been found meet all 
the specifications along the six flight test maneuvers. In the next project phase, these 
controllers will be implemented in a manned simulation of the flight test trajectory control 
system in order to evaluate it more completely. 

The development of a nonlinear controller for a generic aircraft is given in the Ap- 
pendix A. This approach is useful in developing nonlinear controllers for tracking arbitrary 
command histories. 

4.1. Future Work 

The controller development discussed here has a much wider application than the 
flight test trajectory control problem alone. For example, recently there is a strong interest 
at NASA and USAF in cockpit automation. The problem of automatic flight guidance is 
likely to be central issue here. The methods developed in the present work can yield efficient 
schemes for manual or automatic tracking of guidance commands, such as those generated 
by Calise and Moerder [21]. Indeed, one of the early applications of the prelinearization 
technique was in synthesizing an automatic scheme for landing on an aircraft carrier, see 
[22] for details. 

These and other related problems will be of future interest. 
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ABSTRACT 


H 


altitude rate 


Flight test trajectory control systems are 
designed to enable the pilot to follow complex 
trajectories for evaluating an aircraft within its 
known flight envelope and to explore the 
boundaries of its capabilities. Previous design 
approaches were based on linearized aircraft 
models necessitating a large amount of data 
storage along with gain schedules. In this paper, 
the synthesis of nonlinear flight test trajectory 
controllers for a fixed wing aircraft is 
described. This approach uses singular 
perturbations theory and the recently developed 
theory of prelinearizing transforms. These 
controllers do not require gain scheduling for 
satisfactory operation, can be used In arbitrarily 
nonlinear maneuvers, and are mechanized with a 
direct, non-iterative analytic solution. 
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drag force 
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Si » S> St 


Control surface influence 
coefficients for the moments about 
the roll body axis (X) 

Control surface Influence 
coefficients for the moments about 
the pitch body axis (Y) 

Control surface Influence 
coefficients for the moments about 
the yaw body axis (2) 
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altitude 

X body axis moment of inertia 

X-Y body axis product of inertia 

X-Z body axis product of inertia 

Y body axis moment of inertia 

Y“Z body axis product of inertia 

Z body axis moment of inertia 

Component of thrust along the X 
body axis 

Component of thrust along the Y 
body axis 

Component of thrust along the Z 
body axis. 

lift force 

Mach number 

Aircraft mass 

roll body rate 

pitch body rate 

yaw body rate 

total aerodynamic and thrust 
moment about the X body axis not 
Including the control moments 

total aerodynamic and thrust 
moment about the Y body axis, not 
including the control moments 

total aerodynamic and thrust 
moment about the Z body axis, not 
Including the control moments 



5 side force 

T engine thrust 

V total velocity 

a angle of attack 

B angle of sideslip 

6 Aileron Deflection 
a 

6 Elevator Deflection 

e 

6 Rudder Deflection 

r 

# roll attitude 

0 pitch attitude 

• indicates derivative with respect 

to time. 

... an underbar indicates the 

transpose of a vector 

denotes the fast variables in the 
slow-time scale 

denotes the slow variables In the 
fast-time scale 

INTRODUCTION 

The motivation for the development of flight 
test trajectory controllers is well documented in 
the literature 1 . The primary obj ectiv e is to 

enable the ‘pi Tot - to follow complex flight test 
trajectories consistently and accurately. Two 
versions of these controllers have been employed, 
viz a closed-loop automatic system and an open- 
loop system providing manual piloting informatics. 
Originally, the open loop flight test trajectory 
guidance algorithms were developed on-line in a 
piloted simulation using cut and try techniques. 
This approach was not only manpower intensive, but 
often produced less than desirable” controllers. 
Closedkioop system designs based on linearized 
aircraft models ^ required the generation of 
large amounts of numerical data to arrive at 
satisfactory designs. Further, gain scheduling 
was found to be essential for acceptable 
performance . 

The present paper deals with the synthesis 
of nonlinear flight test trajectory controllers 
using the recent results In the pr el I near I zing 
transformation theory! ^ and the singular 
pert ur bat Ions theory 1 . The application of 

singular perturbation theory to this problem 
simplifies the linearizing transformation 
considerably, in addition to providing a 
consistent means for eliminating Ignorable state 
variables. In this frame-work, the sta te 
variables in the original nonlinear problem are 
retained, while the control variables are 
transformed. This Is advantageous from an 
Implementation point of view. Splitting the 
dynamics based on speed of state variable 
evolution generates a controller In which certain 
control loops can be computed at a slower rate 
than others on the flight control computer. It Is 
interesting to note that in Reference 9, even 
though the controller development did not make use 


of singular perturbation theory, the time-scale 
separation formed a basis for implementation on 
the flight control computer. Note, however, that 
the flight test trajectory control problem 
discussed here Is distinct from those described In 
references 7*10, since In those investigations, 
the trajectory to be followed consisted of three 
position components specified as a function of 
time. 

The flight test trajectory controller 
synthesis for a fixed wing high performance 
fighter aircraft without VTOL or hover 
capabilities will be discussed In this paper. 
Though specific engine and airframe models are 
required for implementation, these will not be 
discussed as they are not central to the material 
to be presented. It will be assumed that the 
aircraft under consideration has the four usual 
controls throttle, aileron, rudder and ^levator. 
It is further assumed that the aircraft has no 
direct force generation devices other than the 
engine thrust without thrust vectoring 
capabilities. The task of the flight test 
trajectory controller Is to track the given 
commands In airspeed angle of attack, ^ngle of 
side slip and altitude In the presence of 
disturbances and modeling uncertain! ties 

The next section will discuss the aircraft 
modeling and time-scale separation The details 
of prelinearization and slow-fast controller 
synthesis will be given in Section 3- Simulation 
results using the nonlinear flight test trajectory 
controllers will be presented In Section U, 


MODELING AND TIME -SCALE SEPARATION : 


The equations of motion for an aircraft 
flight over flat, nonrotating earth with zero 
ambient winds is given by 

V • [-D cosB + S sinB ♦ (ki cosa cosB 
♦ k, sina cosB)T -mg (sine cosa cosB 
- cose sin# sinB - cosG cos# 
sina cosB)]/m ( 1 ) 

a - [-L ♦(kiCosa - k l slna)T ♦ mg (cose cos# 
cosa + sine sina)] /mV cosB ♦ q — tanB 
(pcosa ♦ r sina) ( 2\ 


B - [D sinB ♦ S cosB - (k, cosa sinB * k, sina 
sinB)T ♦ mg (sine cosa sinB ♦ cosB sin# 
cosB cose cos# sina sinB)]/Vm 
♦ p sina - r cosa (3) 


h - V(cosB cosa sin6 - sinB sin# cos0 

- cosB sina cos# cosS) (*0 

# - p+q sin# tan© ♦ r cos# tan© (5) 
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9 » q cos# - r sin# 


(6 


where 


ep - fPI, ♦ RI| ♦ e,I, 6 * f* 2 I i 6 + e # I j 6 


♦ pqdjjj 1 ! - D z I.) - qr 

(D x I,4I ** 1,)3/I 


(7) 


eq - [QIj ♦ f,I» « e * f.U « a + f.I. « r 

♦ P !„!* " prD y I.+r I XZ I»]/I 


( 8 ) 


ef - TPIi ♦ RI« + g,I» * e * gi^t * a ♦ g.I. ® p 
♦ Nd I.-OI.) ~qr (D Ii* (9) 

XZ Z X 


wi th 


1,-11 

1 y z 

i, - i i 

y xz * 

i„ - i I -I 

x z xz 

I*-II 
X y 

D - I -I 
x z y 

D - I -I 
y x z 


D - I -I 
z y x 


I - I I I -I I 
x y z y xz 


I -I -0 for the Aircraft under 
xy yz 

consideration 


a 0 - b„ 

a, - b, cos# ♦ b*3in$ tan0 
a 2 - b*, cos# tan© - bjSin# 

and 

b 0 - (cosB cosa sin© - sinB sin# cos© 

-cosB sina cos# cos©) 

b, - V(-sinB cosot sine - cosB sin# cose 

♦ sinB sina cos# cose) 

b 2 - V(-oosB sina Sine - cosB cosa cos# cos©) 

b, - V { cosB cosa cos© ♦ sinB sin# sine 

♦ cosB sina cos# sin0) 

b* - V(-sin0 cos# cos© * cosB sina sin# cos©) 

In order to illustrate the singular perturbation 
procedure, the expressions (1)-(3), {5)-(1l) are 
next expressed in a compact form as 

x - A , (x) + B,(x) z * C, (x) u, Hi, 

ez - A 2 (x,u, , z) ♦ Bj(z) ♦ C 2 (x) u, (13) 


with 


The X, Y positions and yaw attitude # are 
ignorable In the flight test trajectory problem 
under consideration. Consequently the equations 
describing their dynamics hasCbeen eliminated. 

The interpolation parameter c Introduced on the 
left-hand-side of equation (7)-(9) is motivated 
from the forced singular perturbation theory * and 
serves to indicate the difference In time-scale 
between the expressions (l)-(6) and the body rate 
equations (7)- (9). Thus, with e-0, one obtains 
the slow- time scale problem, while e-1 yields the 
complete system. Assuming that the control 
surface deflections have a relatively small effect 
on lift, drag and sldeforce, with e-0, the fast 
variables p, q and r appear ’control like’ in the 
system (1) (6) with three nonlinear algebraic 
equations relating them to control surface 
deflections. 

This approach runs into difficulty however, 
since the expression (H) does not contain p, q, r 
components explicitly. To remedy this situation, 
the expression (*0 Is differentiated once with 
respect to time and substituting for #, $ from 
expressions (5) and (6) one obtains 

h - H (10) 


H - b 0 V ♦ b,B + b 2 a ♦ a 0 p ♦ a,q + a 2 r (1* 


x - [V, a, B, 0 t #, h, H], 


u, - T 


I - [p q r] 


H, - £«. « a « r 3 

A nonlinear controller for the system (12), (13) 
can be designed by transforming it into 
Brunovsky’a canonical form 12 . But this would 
involve the computation of partial derivatives of 
the terms A,(x), B,(x) and C,(x). This difficulty 
is avoided by Invoking the assumption that the 
body rates p q and r evolve faster than other 
variabl es. 

In the following, slow - fast controller 
synthesis using singular perturbation theory is 
discussed without any theoretical development. 
Setting c-0 in the system (12), (13) one obtains 
the slow system as 

x - A , (x) ♦ B 1 (x)z * Cjfx) u, (Hi) 
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0 - A a (x,u lf z) ♦ B a (z) ♦ C a (x) u j (15) 

z are the values of fast state variables in the 
slow-time scale. ^Ideally, as in ref. 16, one 
should solve for z ln/terms of u 2 from the 
expression (15) and substitute in the equations 
(1*0 to obtain a system independent of z. This i* 

difficult in the aircraft trajectory control 
problem due to the^nature of the functions, and 
B a . Alternati vely/nonl inear controller can be 
synthesized for the dynamic system (1ij) with z and 
U| as the controls to track the required x 
commands, see ref. 16 for example. Next the z 
obtained from this^exercise can be substituted in 

(15) to solve for u t provided that C a (x) is 
invertible. This completes the design of slow 
time scale system. 

To derive the fast- time scale controller, one 
assumes that the slow variables are constant in 
the fast-time scale dynamics. Subtracting (15) 
from ( 13 ) and putting 

Az • x-z, Au - U t - U t 

one has 

Ax - A t(x,U| ,Ax) ♦ B > (Ax) ♦ C 2 (x) Au ( 1 6 ) 

x, u, are values of slow state variables in the 
fast- time scale problem, A nonlinear feedback 
controller can again be designed for the system 

(16) to maintain Az close to zero. As long as the 
body rate dynamics remain faster than other 
dynamics, one would expect a satisfactory 
performance from this controller The time scale 
separation similar to the one discussed here has 
been employed in the past for flight control 
system design and is known to be valid in most 
situations. If the actuator dynamics are to be 
included in the control system synthesis, they can 
be handled in an additional time scale. Thus, 
with time scale separation, the flight test 

* ■ i. ectory controller will be of the form given in 
r re 1 , 

In summary, the singular perturbation schen> 
described relies on the fact that the force 
generation mechanisms on the airframe are slower 
than moment generation processes, 

NONLINEAR CONTROLLER DESIGN 


The main specification on the flight test 
trajectory controller is that it should track 
given time histories of V, h and a. While there 
are no direct specifications on the angle of side 
slip 0, it is desirable to maintain it close to 
zero throughout a given maneuver. In the 
following the slow and fast time scale 
controllers for tracking these will be discussed 
separately with specific details. 

slow- time Scale Controller: 

Setting the interpolation parameter to zero 
the slow-time scale dynamics given by the 
expressions (1)-(3), (5), (6) and (10), (11) can 
be written as 

V - c 0 * c,T 


c 0 - [-DcosB ♦ S sinot - mg (sln0 cosa cosB 
- cosG sin$ sinB - cos© cos$ sina 
cos$)]/m 

Cj - (kjcosa cosB ♦ k,sina cosB)/m 
and 

o • c, + c *T ♦ c*p ♦ c q ♦ c*r (18' 

with 

c a - [-L ♦ mg (cos© cos$ cosa * sine 
sina)]/nfifcosB 

c t - (k, cosa - k,sina)/iBVcosB 
c„ - -cosa tanB 
c* - 1 

c 6 - -tanB sina 
and 

6 - d 0 ♦ d,T +d A p +d,r (19) 

with 

d 0 - CD sinB + S cosB + »g (sine cosa sinB 
+ cos6 sin$ cosB - cos© cos$ sina 
sin6)]/Vm 

dj - (-k, cosa sinB - k t sina sin0)/Vm 
d 2 - sina 
d s » -cosa 

♦ - * d*q + d*r (20) 

§ - d*q ♦ d,r (21 ) 

w^h 

d„ * 1 

d» - sin$ tan© 

d« - cos $ tanB 


(17) - cos* 

d, • - 
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‘Sin* 
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fi - H (22’, 


» •* ( • • *» 

H - h • b 0 V ♦ b|B + b,o * a e p 

+ a,q *a 2 r (23) 

0 ■ [PIj * RX « + e t Ii 6 + o 2 X ) 6 * e,I , 6 

^ v o r 

♦ P qd X2 I i ' Dj Ii) ■ g r 
(D x Il+ 1 xz 1,)]/I (211 


o - [qi* ♦ f,i. ? e ♦ r t i„ « a ♦ f,i. « r 

♦ p*I X2 I, - p rD y I. ♦ r*l n UVl (25: 

0 - [PI, ♦ RI, * 8,1 ( g,I* « B * g.I, i r 

* P qd X2 I»"I> 8 I«) -q r 

(DI,‘ I„I.)]/I (26) 

X xz 


The system (17)“ (26 ) describes the slow-tlrae scale 
dynamics with Control variables p, q, r and T, the 
body rates and the engine thrust. Since there are 
seven state variables and four control variables, 
only four of these states can be completely 
controlled. 

The variables of interest in the present 
flight test trajectory control problem can be 
broadly grouped into two sets i.e., tracking the 
variables 


It can be verified that except for the altitude 
loop, all other controllers will have zero steady 
state errors for ramp commands. Moreover, the 
natural fre quencfc* and damping ratios of these 
control loops are 

(i) airspeed loop: 

k, 

“n v - ^ • *V ' VCf 

(ii) Angle of attack loop: 


k. 



(iii) Angle of side slip: 


k* 



(iv) Altitude loop: 

k* 

“a ‘ / * r «h ’ 2/ilf 


Hence, two time response specifications can be set 
for each of these control loops. Note, however, 
that the extent to which these specifications are 
met would depend on the actuator saturation levels 
and the tracking commands. 


(i) 

v, h, a, 6 

(ii) 

V, h. 6. * 


Consequently, the number of states to be 
controlled is equal to the number of controls. 

The other state variables are treated as free. 

The discussions in the following will be limited 
to Case (1). Case (ii) can be handled In an 
entirely analogous manner. Since there are four 
i tat e- variables to be tracked and four control 
variables, the slow-time scale system can be put 
in the following form by defining four new pseud 
controls U,, U a , U t , U*. 

h - H 

H - U, 

V - U* (27) 

a - U, 
i - U„ 

Next, four independent linear controllers can be 
designed for this system to ensure zero tracking 
errors for ramp commands. Typically, the 
airspeed angle of attack and angle of side slip 
will have proportional plus Integral control while 
the altitude loop will have a proportional plus 
derivative control. The form of these pseudo 
controller loops amm^given in Figure 2 and 3. 


The remaining task in the slow-time scale 
control problem Is that of converting the four 
pseudo controls to the four real controls. Thus, 
from (17)-(19) and ( 23 ), one has 


T - (U a - e.)/c, 


( 28 ) 


a O a 2 


P 


^-(boU, ♦ b,U„ * b 2 U, ) 

c„ c* c* 


Q 

m 

U 2 c, _ CjT 

d 2 o d| 


r 


U„ d 0 - dj 

» 


. . 


p m 


(29) 


It is assumed here that the required thrust givrr 
fcy the expression (28) can be converted into the 
throttle setting using tabular data. The set of 
linear equations (29) can be solved for p, q, r. 

If the aircraft Is in symmetric flight, the 3x3 
matrix premulti plying p, q, r vector will have 
less than three rank. This is because under these 
conditions, In the reduced order problem, there 
are three longitudinal state variables to be 
controlled while there are only two available 
controls, viz; the thrust T and the pitch body 
rate q. Hence one can precisely control the 
airspeed and altitude or the airspeed and angle of 
attack or the angle of attack and altitude. 
However, if a linear combination of any two of 
these three variables amm/formed, then that linear 
combination can be controlled exactly along with 
the remaining variable. The situation In the 
lateral channel Is reverse under the symmetric 
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flight condition, l.e., there are two control 
variables, r and p available while there is just 
one state, B, to be controlled. Thus, in the 
lateral channel, the angle of side slip can be 
maintained zero either by the yaw body rate r or 
the roll body rate ? or by a linear combination of 
these. The approach adopted In the present work 
is to use r to control B while using p to maintain 
the roll attitude zero. 

Once the expressions (28) and (29) have been 
solved, the p, q, r values and thrust can be 
substituted in equations ( 2 ^-( 2 £) to compute the 
control surface deflections $ r along the 

"outer" solution. • 

Fast- time scale controller: 

In the fast-time scale, the slow-time scale 
variables h, V, a, B are assumed to remain 
constant. The fast-time scale controller attempts 
to maintain the body rates p, q, r close to their 
values in the "outer" solution. Subtracting the 
expressions (2*0- (26) from ( 7 )-( 9 ), and putting 

tm mm tm 

Ap - p-p, Aq - q-q, Ar - r-r 


one obtains 


there are no explicit tracking requirements. The 
real controls can be obtained from pseudo controls 
as 


e,I, e,I, e,I, 


(4 e - V 


rS 

f.I. f.l. M. 


(S a - V 

m 

fs 

81^* Bi^« 8»J * 


(6 r - V 

■ 

r 

* 8 

• « 




• 


(3*0 


where 

f„ * IU, - (P P)I, - (R-R)I, - (Ap • Aq <► p*Aq 

♦ q-Ap) (I^I, - D Z X «^ " (Aq -Ar 

♦ q • Ar ♦ r-Aq) (D I, + I T s ) 

H H x 1 xz 

f* - IU § - (Q - Q)I* - (Ap ♦Ap p) I xz I, 

+ (Ap*Ar + Ar-p ♦ Ap-r) D I„ 

2 / 

- (Ar ♦ Ar«r)I I* 

f* - IUt (P-P)I, - (ft-R)I* ♦ (Aq Ar ♦ q-Ar 

♦ r*Aq)(D x I,+ 


Ap - [(P-P)^ ♦ (R-R)I, ♦ e,I,( 6 e - 6 ) 

♦ e 1 I 1 ( 6 a - S a ) ♦ e,I , ( 6 p - 6 r > ♦ 

(Ap-Aq ♦ p-Aq + q Ap) (I^Ij - D I,) 

- (Aq* Ar + q Ar + r-Aq) (D^I, ♦ 

(30) 

Aq - [CQ-Q)I % 4 f 1 I <r ( 6 e - 5^) ♦ - 6 ) 

+ “ « r ) * (Ap 4 Ap-p) I^I* 

- (Ap*Ar ♦ Ar*p ♦ Ap*r) D y I* 

♦ (Ar 4 Ar-r ) ( 31 ) 

Ar - f(P-P)I, 4 (R-R)i, ♦ g 1 I t ( 6 e - l e ) 

♦ g*I«(A a “ 6 ft ) 4 gl i t ( 6 r - A p )- (Ap*Ar 

♦ q-Ar 4 r-Aq) (D^I t + I xz I # )3/I (32) 

The fast- time scale controller maintains Ap, Aq, Ar 
close to zero throughout a given maneuver* Since 
there are three Independent controls available, 
three pseudo controls are next defined such that 
the system C 30 )- C 32 ) takes the form 

Ap * U, 

Aq - U 4 (33 

Ar - U 7 

Three independent control loops can be designed in 
pseudo controls U f , U 4 , U? such that the system 
(33) has a much faster time constant than the slow- 
time scale system (17)- (26). Integral feedbacks 
are not necessary for these control loops since 


The set of linear algebraic equations (3*0 can be 
solved for (6 -6 ) , (A -A.) (A ” 6 ^). Note that 
the matrix muftiplyingthese quantities has full 
rank everywhere on the flight envelope and a unique 
solution always exists. Since A , A , 6 are known 
from the outer solution, the actual control su^fac- 

0<r ' r • * - nr r -iir t " 

NONLINEAR CONTROLLER EVALUATION 

In order to test the performance of the 
nonlinear controller synthesized in Section 3, it 
is implemented on a slx-degrees-of-freedora 
simulation of a higft performance fighter aircraft 
including a first order engine dynamics. Two 
symmetric flight test maneuvers were executed, viz, 
a level acceleration trajectory and a pushover- 
pullup trajectory. 

The level acceleration trajectory Is a wings 
level, constant altitude maneuver with a ramp Mach 
number command. This maneuver is Initiated at 
20000' and Mach 0.9. The objective Is to 2 
accelerate the aircraft at about 1 0 ft/sec for 5 
seconds. An initial constant Mach number leg is 
included in the command history to provide adequate 
time for damping out the effect of inexact trim 
conditions. The Mach number command as well as the 
response of the aircraft are given in Fig. The 
altitude history is given in Fig. 5, and the 
throttle setting is in Fig. 6 . The Mach nunber 
tracking error is less than ^ 0.001 during most of 
the acceleration phase. The altitude is maintained 
within ^0.2 feet. This maneuver required the 
after-burner thrust as can be seen from the 
throttle history in Fig. 6 . 

A pushover- pull up trajectory Is executed next. 
This flight test trajectory is a wings level, 
constant Mach nunber maneuver in which the angle of 
attack is varied a specified Increment about the 
trim value at some specified rate. Fig. 7 depicts 
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the commanded and the actual angle of attack. Th< 
angle of attack tracking error Is well within 
jj;0.1 o , The Mach number was maintained within 
♦0.006 throughout the maneuver as can be observed 
from Fig. 8. A transient appearing at 20 seconds 
in the Mach number history is due to the throttle- 
thrust characteristics of the engine. For this 
aircraft engine, there Is small core thrust 
saturation region approximately between 83 and 98 
degrees of throttle setting. This region is 
illustrated on the throttle setting curve in Fig. 

9. Note that in this maneuver, the altitude is not 
controlled and it is free to change as it may 
throughout the trajectory. In Fig. 10, the 
resulting altitude is given. During the initial 
negative angle of attack rate region as well as for 
some positive angle of attack rate, the aircraft 
loses as much as 2000 feet. During the second half 
of this maneuver, the aircraft gains altitude and 
overshoots the Initial altitude by almost 700 feet. 
If desired, the commanded angle of attack history 
can be tailored to return the aircraft to initial 
straight and level conditions at the end of the 
maneuver. 

Summarizing the results presented so far, the 
nonlinear controller performance for these two 
maneuvers has been found very good. Currently, 
work Is underway to evaluate the controller 
perf orifice along several other maneuvers. 


CONCLUSIONS : 

Nonlinear controllers for tracking flight te> * 
trajectory commands were described. The controller 
development employed singular perturbation theory 
a ~nd the recent results for a class of nonlinear 
systems. The synthesized controllers do not 
require gain scheduling and can be implemented v 

different rates on the flight control computer. 

The approach presented here can be extended for 
tracking three positions components and velocity 
for example, as in an autoland task. The 
synthesized control laws have a general -haracter 
in the sense that their form remains invariant for 
any conventional aircraft. However, the time-scale 
separation presented must be valid for the maneuver 
under consideration. 
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Figure A. Mach Number Evolution Alon^ The 
Level Acceleration Flight lest 
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Figure 5. Altitude History Along The 
Level Acceleration Flight 
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Figure 6. Throttle Setting Along The 
Level Acceleration Flight 
Test Trajectory 

A: Core Thrust 
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Figure 7. Angle of Attack Along 
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Figure 8. Mach Number Evolution 
Along The Pushover- 
Pullup Flight Test 
Trajectory 
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Figure 9. Throttle Setting Along The 
Pushover-Pul lup Flight Test 
Trajectory 
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Figure 10. Altitude Evolution Along 

The Pu ah ov e r-Full up Flight 
Test Trajectory (seconds, 
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